Riesz-type inequalities and maximum flux exchange flow 



I McGillivray 
School of Mathematics 
University of Bristol 
University Walk 
Bristol BS8 1TW 
United Kingdom 
e maiemg@bristol.ac.uk 
t + 44 (0)117 3311663 
f + 44 (0)117 9287999 



Abstract 



o 

> . 

Q , Let D stand for the open unit disc in R d (d > 1) and (D, S3, m) for the usual Lebesgue measure 

■ space on D. Let J4f stand for the real Hilbert space L 2 (D, m) with standard inner product 

(■, ■). The letter G signifies the Green operator for the (non- negative) Dirichlet Laplacian 
—A in .jtff and ip the torsion function G\d- We pose the following problem. Determine the 
optimisers for the shape optimisation problem 

Pi , . at := sup \ (Gxa,xa) ■ A C D is open and (ip, xa) < t\ 

. , where the parameter t lies in the range < t < (ip, 1). We answer this question in the one- 

dimensional case d = 1. We apply this to a problem connected to maximum flux exchange 
flow in a vertical duct. We also show existence of optimisers for a relaxed version of the above 
variational problem and derive some symmetry properties of the solutions. 
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1 Introduction 

CO 

Let S! stand for a bounded open set in M d (d > 1) and (CI, £§, m) for the usual Lebesgue measure 
space on CI. Let J$? stand for the real Hilbert space L 2 (Cl, m) with standard inner product (•, •). 
The letter G signifies the Green operator for the (non- negative) Dirichlet Laplacian — A in 34? and 
ip the torsion function Gxn- We pose the following problem. Determine the optimisers for the 
shape optimisation problem 

X 1 a t := sup { (Gxa,Xa) ■ A C Q is open and (ip,XA) < t \ (1.1) 

where the parameter t lies in the range < t < (ip, 1). We show that optimisers exist for a relaxed 
version of this problem and derive certain symmetry properties of the solutions when Q is replaced 
by the open unit ball D. We obtain the explicit form of the optimisers in the one-dimensional 
case d — 1 for the open interval D — (— 1, 1) . 

Define 

V t := | / e3T : < / < 1 m-a.e. on CI and (/, ip) < t \ 
for t in the range < t < (ip, 1) and consider the relaxed variational problem 

ft:=su P { J(f): feVt}, (1.2) 
where </(/) = (/, G/). The first main result runs as follows. 
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Theorem 1.1. For each t in the range < t < (ip, 1), there exists f £ Vt such that f3 t = </(/)• 



In case f2 is replaced by the open unit ball D centred at the origin, we can say more about the 
symmetry properties of optimisers. In fact, 

Theorem 1.2. Let f £Vt such that /3 t = </(/). Then f possesses circular cap symmetry. 

We now turn to the one-dimensional case d — 1 so that D = (—1,1) and the torsion function 
tp is given explicitly by ip(x) — (1/2) (1 — x 2 ) for x £ D. Noting that (ip, 1) = 2/3, define 
tp : D — > (0, 2/3) by ip(x) :— (xtx,i)i V0j an d specify £t £ (—1, 1) uniquely via the relation 



for each t £ (0, 2/3). Set A t := (&, 1). Then 

Theorem 1.3. For any open subset A in D satisfying (ip, xa) < t it holds that 

(Gxa, Xa) < (GxA t , XA t ), 
and equality occurs precisely when either A = A t or A = —A t . 

The inequality is somewhat reminiscent of the Riesz rearrangement inequality: this justifies the 
epithet in the title. This problem has a probabilistic interpretation in so far as the function G xa 
is the expected occupation time in A spent by absorbing Brownian motion in D (associated to 
the Laplacian A). The d > 2 case has not yet been resolved. It is tempting to speculate that a 
hyperbolic cap optimises (jl.ip in this case. Numerical evidence does not seem to bear this out, 
however [S]. 

One reason why this problem is intriguing is because of its connection to maximum flux exchange 
flow in a vertical duct, a model of lava flow in a volcanic vent (see [4]). In the two-dimensional 
case d = 2, we imagine a configuration of two immiscible fluids in D x R with different physical 
characteristics in a state of steady flow. The densities of the fluids are labelled p, p' and we take 
p > p' . Each fluid has unit viscosity. With respect to cylindrical coordinates (x, z) £ D x R, 
gravity acts in the direction (0, —1) according to the model. The pressure p depends only upon z 
and has constant gradient dp/d z — —G. Suppose that the fluid with density p occupies a region 
in D x R with cross-section A C D. Restricting the problem to D, the velocity u of the components 
of the fluid may be described (informally) using the Navier-Stokes equation via 

= Au + G - pg on A; 
= Au + G-p'g on D\A. 

Non-slip (Dirichlct) boundary conditions are imposed on the boundary of D. It is also assumed 
that u and its gradient are continuous on the interface between the two regions A and D \ A 
(continuity of velocity and stress). 

The parameter G lies in the interval (p' g, pg). This allows the possibility of a bi-directional flow. 
Upon rescaling (and relabelling the velocities) we obtain the system 



(1.3) 








= A u — A — 1 on 
= Am-A + 1 on 



A; 

D\A: 



(1.4) 



where 



A : 



(p' + p)g-2G 
(P ~ P') 9 



e(-l 



i) 
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is a proxy for the pressure gradient. Two problems arise. One is to maximise the flux Q := 
(Xd\a, u ) amongst all regions A which satisfy the flux balance condition (u, 1) = with constant 
A; the other in which we optimize also over A. In detail, we seek optimisers for the problems 

7 := sup{ {Xd\a, u) : (u, 1) = 0, A C D open, A G (-1, 1) } , (1.5) 
7a := sup{ (xd\a, u) : (u, 1) = 0, A C D open} , (1.6) 

where in the latter A is fixed in the interval (—1,1). It turns out that problem is closely 

related to the two problems above. Note too that these problems have obvious analogues for the 
case d = 1 . 

We come to our last main result. Note that the d = 2 analogue is discussed as a marginal case in 
®- 

Theorem 1.4. In case d = 1, 

(i) for each A G (— 1, 1), the problem il.6}) is optimised precisely when either A — A i-x or 
A = -Ai-x ; 
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(ii) the problem U.5]) is optimised precisely when either A = (0, 1) or A = (—1, 0) and has 
optimal value 1/12. 

We give a brief sketch of the organisation of the paper. In Section 2, we obtain existence of 
optimisers for the relaxed problem (|1.2j) and derive some symmetry properties when fl is replaced 
by the ball D. Sections 3 to 8 deal with the proof of Theorem 1.3. Section 9 contains an application 
to maximum flux exchange flow (Theorem 1 1.4[) . 

2 Existence of optimisers and symmetry in a general relaxed setting 

For the sake of clarity, we first of all remark that the (non-negative) Dirichlet Laplacian (D(— A), —A) 
is associated with the Dirichlet form (J^, S) in .J? with form domain & := W ' (ft) and 

S (u, v) = (Du, Dv) (u,v€&). 
We begin with the proof of Theorem 11.11 

Proof of Theorem ll.ll Let {f n )neN be a maximising sequence for /3 t . Now, V* is weakly sequentially 
compact in Jf. This follows by appeal to Theorem 10.2.9 due to the fact that V* is bounded, 
closed and convex in the reflexive Banach space Jff. So we may assume that (/„) converges weakly 
in Jtf? to some / G Vt as n oo after choosing a subsequence if necessary. 

Put u n := Gf n . Then for each n, 

IKIIw^cn) ^ \/2m(ft) IM| w i,2 (n) . 
Additionally, 

IWIwo 1 ' 2 ^) = S ( U n, U n ) + (u n , U n ) = (fn, Gf n ) + (Gf n , Gf n ) < (l,1p) + (lp,%l>). 

In short, the sequence (u n ) is bounded in W ' (f2). In case d > 2 by the Rellich-Kondrachov 
compactness theorem ([3] 5.7 for example), we may assume that (u n ) converges in L 1 (J7, m) to 
some element u after extracting a subsequence if necessary. In case d = 1, we use Morrey's 
inequality (see [5] 5.6.2, for example) and the Arzela-Ascoli compactness criterion to extract a 
uniformly convergent subsequence. The details are described in the proof of Theorem 13. II 

For each neN, 

(u n , ip) = {Gf n , if) = (f n , Gf) for all f eJf, 
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which yields 

(u, <p) = (/, G<p) = (Gf, if) for all <p e & 
upon taking limits. Therefore, u — Gf m-a.e. on tt. Moreover, 

J(f)~ J (fn) = {u, f)-(u n , f n ) = (u, f - f n ) + (fn,U-U n ) 

and the right-hand side converges to zero as n — > oo in virtue of the weak respectively strong 
L (CI, m) (or uniform in the case d = 1) convergence of the sequences (/„) respectively (u n ). As 
j3 t = lim ?woo J(f n ) it follows that f3 t = J(f). □ 

In the remainder of this section, we replace CI by the open unit ball D in M. d centred at the origin. 
We first discuss the operation of polarisation for integrable functions on D (see [5] and references 
therein). For v e S^ 1 the closed half-space H = H v is defined by 

H v := {x 6 R d : x-v> 0} 

with an associated reflection 

t h : R d -> R d ; x H- x - 2 (x • v) v. 

Refer to the collection of all these closed half-spaces by %. The polarisation fjj of / € L+(D, m) 
with respect to H e H is defined as follows. Choose an m-version of /, which we again denote by 
/. Set 

f ( \ -/ f( x ) A f( T H x ) for xeDHH, 
lH[X) ' 1 /(x)V/(t h z) for xGi?\ff. 

Its m-equivalence class is the polarisation of /. The definition is well-posed. 
The Green kernel G(x, y) is given by 

G(x, y) - - x) - ^(IxKy - x*)) for (x, y) E D x D \ d, 

where $ is the fundamental solution of Laplace's equation in M. d , d stands for the diagonal in 
D x D and the decoration * refers to inversion in the unit sphere. We note the inequality 

G(x, y) > G(x, r H y) for any x, y e D n int H, (2.1) 

which follows from the strong maximum principle. 

Theorem 2.1. Let f £ L+(D, m) and H E H. Then J(f) < J{Jh) with equality if and only if 
either f = fn or f o th — fn m-a.e. on D. 

Proof. We work with an m-version of /, again denoted /. Define 

A+ :={xeDnH: f(x) < f{ TH x)} 

and similarly B + but with the strict inequality replaced by the sign >. Put A~ := thA + and 
A := A+ U A~. Set S := D \ A. In this notation, 

iff = XA.f°T H + Xs /■ 
As a consequence, 

J(f H ) = J(XA f°r„) + 2 (xa fo-TH, G X sf) + J(xs f) = J(XA f)+2 (xa f°T„, G X sf) + J(xs f) 
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and a similar identity holds for J(f) but without composition with reflection. We may then write 
that 

J(fa)-J(f) = 2 f f (f(THx)-f(x))(g(x,y)-g(T H x,y))(f(y)-f(r H y))m(dy)m(dx). 

JA+ JB+ 

It is clear from this representation with the help of (|2.1[) that J(f) < J(fn)- 

In the case of equality, it holds that either m(A + ) = or m(B + ) = 0. In the former case, f = Jh 
while in the latter, / o rjj = fu m-a.e. on D. 

The spherical cap symmetrisation (see [7], [5], [5] for example) of A £ 38 with respect to the 
direction to G S^ 1 is the set A* G 38 specified uniquely by the conditions 

A*n{0} = An{0}, 

A*ndB(0,r) = B(ru, p)ndB(0, r) for some p > 0, 

oy(A* n <9B(0, r)) = a r (B(ruj, p) n 95(0, r)), 

for each r G (0, 1). Here, oy stands for the surface area measure on dB(0, r). The spherical 
cap symmetrisation of / G ^+(5, wi) (denoted /* for brevity) is defined as follows. Choose an 
m- version of /, which we again denote by /. Let /* be the unique function such that 

{/* >t} = {f> t}* for each t S K. 

Its m-equivalence class is the polarisation of /. The definition is again well-posed. We also write 
/* as C u f. 

Before proving Theorem 11.21 we prepare a number of lemmas. We first discuss a useful two-point 
inequality. We introduce the notation 

Q := {(aii, x 2 ) e M 2 : Xl > and x 2 > 0} , 
R := {(xi, x 2 ) e Q : < x 2 < x\} , 
S := {(xi, x 2 ) G Q : < xi < x 2 } . 

Equip Q with the ^-norm ||a;||i := \xi\ + \x 2 \ where x = (xi, x 2 ) G Q. Define a mapping 
tp : Q — >■ Q via 

(xi, x 2 ) H 1 (xi V x 2 , xi A a; 2 )- 

A geometric argument establishes the following lemma. 

Lemma 2.1. For any x, y G Q, 

Wfx - ipy\ x < \\x - y\\i 

with strict inequality if and only if x € R and y G S or x G R and y G S or the same with the 
roles of x and y interchanged. 

For lu G S f<i ~ 1 introduce the collection of closed half-spaces 

H u := {x G R d : x ■ v > 0} . 
Lemma 2.2. Let / G L\(D, m) and lu G S d 1 . For any H G % w , 

WfH - Cw/||ii(£),m) < 11/ - Cw/Hlha™) ( 2 - 2 ) 
wii/i strict inequality if 

m({f ot h > /}) > 0. 
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Proof. Select an m- version of /, again denoted /. Note that = /*. By the two-point inequality 
Lemma 12.11 

\f H (x) - f*(x)\ + \f H (r H x) - f*{r H x)\ < \f(x) - f*(x)\ + \f{r H x) - f*(r H x)\ (2.3) 
for x G D n H . It only remains to integrate over D n H to obtain the inequality. 

For each x G D n H the pair (f*(x), J*(thx)) belongs to R. By Lemma 12.11 the condition 
(f*(x), J*(thx)) G S guarantees strict inequality in (|2.3[) . This observation leads to the criterion 
in the Lemma. □ 

The next lemma is a spherical cap symmetrisation counterpart to [2] Lemma 6.3, and extends [7] 
Lemma 3.9. 

Lemma 2.3. Let f G L+(D, m) and uj G S d ~ x and assume that f =^ C u f. Then there exists 
H G T~L U such that 

\\/h - C u f\\ L i(D,m) < 11/ _ Cw/IUi^.m)- 

Proof. For shortness, write /* for C u /. As / ^ /* there exists t > such that 

m({f>t}A{f*>t})>0. 

It follows that the sets A := {/ < t < /*} and B :— {/* < i < /} are disjoint and have identical 
positive m-measure. 

We claim that there exists H G such that m(A n T//i3) > 0. Taking this as read, on A n TgB 
we have that /* > t > f * o th so that A n TgB C H . Also, / < t < f o th there. In short, 
A n t h B C{/or ff > /} n H. So m({f o r H > /}) > and there is strict inequality in (j2"T2"|) by 
Lemma 12.21 

To prove the claim, we assume for a contradiction that m(A n thB) = for all H G Let f 
be a countable dense subset in S d ^ 1 n Then 

m(in (J r^B) = 0. 

Therefore, for all r G (0, 1), it holds that 

a r (A r n TH„B r ) = for every v G F, 

except on a A- null set N. Here, A stands for Lebesgue measure on the Borel sets in R, and 
A r := A n 8B(0, r) for the section of A (likewise for B r ). Let ^ G S' d_1 PI i? w with corresponding 
reflection r = th v - Select a sequence (fj) in F which converges to v in S d ^ 1 . Write Tj for the 
reflection associated to closed half-space H Vj . For r G (0, 1) \ N, 

\a r (A r r\rB r ) - cr r (A r H TjB r )\ < \\xb - Xb ° r o T,-|| L i(a B ( ,,.),a r ), 

and this latter converges to zero as j — > oo. This is due to the fact that the special orthogonal 
group SO(d) acts continuously on L 1 (S d ^ 1 : a). We derive therefore that 

a r (A r n T H „B r ) = for every v G S^ 1 n (2.4) 

for all r G (0, 1) \ N. 

To conclude the argument, choose r G (0, 1) \ N such that a r (A r ) = a r (B r ) > 0. Use Lebesgue's 
density theorem to select a density point x for A r lying in A r , and choose y in B r similarly. Then 
f*(x) > t > f*(y). So there exists v G S d ~ x n such that with r = th k we have that ry = x. 
But this means that 

,. a r (A r DtB,. n B(x, e)) 

lim = 1 

ao a r (9B(0, r)f!5(i, e)) 
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so that, in fact, a r (A r P\rB r ) > 0, contradicting (12.41) . □ 

Proof of Theorem \1.2\ Assume for a contradiction that / ^ C u f for each ui 6 £ Then there 
exists to G such that 

<5 := in| ||/ - C v f\\ L xm rn) = \\f - C u f\\ L i/ Dm ) > 0. 
By Lemma \2. 31 there exists H G H w such that 

WfH - C u f\\Li(D,m) < 11/ _ C u f\\LHD,m)- 

It is plain that f fn- But also / o tjj ^ f 'h , for otherwise, 

11/ - C auJ f\\ L i( D:m ) = \\fu - C^/Hl^D,™) < 11/ - Cw/|Ui(£),m)) 

contradicting optimality of w. It follows by Theorem 12.11 that J(f) < J(fii) and this contradicts 
the optimality of / in the expression for f3 t . □ 

3 Preliminaries for the one-dimensional problem 

In the remainder of the article we work in the one-dimensional setting where D = (—1, 1). In this 
context, the corresponding Green operator G has kernel given by 

for x, y G D. We record the useful inequality 

| G(x, y) — G(x, x)\ < | y — x | for all x, y G D, (3.2) 

for future use. As noted above, the torsion function ip := G\d is given explicitly by ip(x) = 
(1/2) (1-x 2 ) forxGA and 

(l,V0 = 2/3. (3.3) 
The Green kernel may be bounded in terms of that is, 

G(x, y) < i}){x) for all y G D, (3.4) 
with fixed 

For t G (0, 2/3) introduce the shape space 

U t := { / = XA ■■ A C D is open and (/, V) < *}• 
We may then write 

at = sup{ J(f) : /G t/ t }. (3.5) 

For each i G (0, 2/3) and m G N define l/ t to be the collection of all functions of the form / = \A 
where A is a union of at most m disjoint open intervals in D with the additional requirement that 
{f, "0) < t- We occassionally refer to the condition 

intA = A. (3.6) 
We also introduce the variational problem 

a< m) :=sup{j(/): / € !7 t (m) }. (3.7) 
We now derive the crucial property that p. 71) attains its optimum. 
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Theorem 3.1. For each t £ (0, 2/3) and m 6 N t/iere exists f G uj: m ' with (/, ip) — t such that 

* { r } = j(f). 

Proof. Let (/ n ) n gN be a maximising sequence for a[ m ' ) . Each /„ may be written in the form 
fn = Z)j=i XA nj for some 1 < k n < m where A„j = (a nj , b nj ) and 

-1 < a n i < b n i < a n2 < b„2 < • • • < a„t„ < 6„ fcn < 1. 

After selecting a subsequence if necessary we may suppose that k n takes a fixed value k for some k 
between 1 and m. On appeal to the Bolzano- Weierstrass theorem, we may assume (perhaps after 
discarding a subsequence) that a n j — > dj and b n j — > bj as n — >• oo where 



Set / := Ylj—i XAj where Aj — (aj, bj). By the dominated convergence theorem, (/„) converges 



-1 < oi < bx < a 2 < b 2 < • ■ ■ < a k < b k < 1. (3.8) 

Ek 
j = l X- A 3 

weakly to / in 3? . 

Put u n := G/„ as before. Then the sequence (u n ) is bounded in W ' (D) as in the proof of 
Theorem ll.il By Morrey's inequality (see [3] 5.6.2, for example), 

|| u n \\ c o,i/2m) < c f° r all n G N 

for some finite constant c; in particular, 

| u n (x) - u n {y) | < c | x - y | 1/2 for any x,y € D 

and any neE Thus, (u„) forms a bounded and equicontinuous sequence in C(D). By the Arzela- 
Ascoli compactness criterion, we may assume that (u n ) converges uniformly to some u G C(D) as 
n — > oo after extracting a subsequence if necessary. Now continue the argument as in the proof of 
Theorem 1 1.1 1 to conclude that ctj™^ — J(f). 

We now show that (/, ip) = t. First note that (/, ip) < t; this flows from the fact that / is a 
weak limit of elements in £/ t . Suppose for a contradiction that (/, ip) < t. As (/, ip) < 2/3, 
in (|3.8[) there must exist j = 0, . . . , k such that bj < Oj+i with the understanding that &o := ~ 1 
and a k +i '■= 1. By choosing B to be a suitable (semi-)open interval in [bj, Oj+i] we can arrange 
that the function fi := f + \b satisfies the requirement (/i, -0) < t as well as J(f) < J(fi). This 
contradicts the optimality of /. □ 

We now revisit the operation of polarisation in the one-dimensional setting. We use the letter P 
to signify the polarisation operator with respect to the closed half-space [0, oo). Thus, for / £ Ut, 
the polarisation is defined by 

Pf(x)-{ fW^fi-x) if 0<x<l, 
I[ ' '~ 1 f(x) Af (-x) if -Kx<0. [6 ^> 

Alternatively, suppose that / = \A where A is an open subset of D. Then Pf = xpa where PA 
denotes the polarisation of the set A; in other words, 

PA = ADtA\J(ALItA) n(0, 1) (3.10) 

where r : D —5- D stands for the reflection x n- —x. We shall sometimes refer to the symmetric 
resp. non-symmetric parts of PA; that is, 

A x := ADtA; 

A 2 := (A U tA) n (0, 1) \ A n tA. l ° > 



Lemma 3.1. Let f G Ut for some t G (0, 2/3). Then the following statements are equivalent: 
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(i) f e PU t ; 

(ii) f = lonS:={xG (0, 1) : f{-x) = 1 }. 

Proof. Let / £ Pi7 t so that f = Pg for some g € C/ t . Let x G (0, 1) with f(-x) = 1. Then 
1 = /(-x) = Pg(-x) = g(x) A ff(-x). So g(a;) = 1 and f(x) = Pg{x) = g(x) V g(-x) = 1. On 
the other hand, suppose that / = 1 on S. For x 6 5, 

P/(x) = 1 V /(-or) = 1 = /(a) while P/(-x) - 1 A f(-x) = f(-x), 

and for x G (0, 1) \ 5, 

P/Or) = /(*) V = f(x) while P/(-x) = /(*) A - - f(-x). 

In other words, / = Pf. □ 

It is sometimes useful to polarise with respect to the closed half-space (— oo, 0]. To distinguish 
between these two polarisations we use the notations P+ , P_. In particular, 

P fW ._f /(^A/H if < a; < 1, , , 

^- J[X) - 1 f{x) V/(-x) if -l<a:<0. 1 j 



Lemma 3.2. Let f = \A £ P+£7 f /or some to G N and t G (0, 2/3) where A satisfies condition 
A3. 6]) . Put g :— \B where B :— D \ A. Then g is an m-version of 1 — / and g G P_C/ 3 / 2 -t- 

Proof. We may suppose that A = Uj=i Aj f° r some 1 < k < to and A, = (oj, fej) with 

— 1 < ax < b\ < a2 < 62 < • • • < a*; < bf~ < 1. 

We use the criterion in Lemma f3.ll Let x G (— 1, 0) such that g{— x) = 1. We first note that a; 
cannot be a boundary point (that is, x $ {di, . . . , afc, 61, . . . , bf. }). For if it is, then cither — x is a 
boundary point or —x £ A. This is due to the fact that / is polarised to the right. In either case, 
we obtain the contradiction that g{— x) = 0. We want to show that g(x) = 1 so suppose on the 
contrary that g(x) = 0. Then for y = — x £ (0, 1), it holds that f{—y) = 1, but f(y) = 0. This 
counters the fact that / £ P+Ut by the criterion. □ 



4 A (non-)optimality criterion 

In this section we develop a (non-)optimality criterion for configurations / in J7 t . Given / G Ut 
define u := Gf. It is known that D(-A) = Wq' 2 (D) n W X ' 2 {D). Thus, u G iy 2 ' 2 (_D) and by a 
Sobolev inequality (see [3] 5.6.3 for example), u belongs to the Holder space C 1 ' 1 / 2 ^). Define 




Then h G C(D) and by l'Hopital's rule, 

h(-l)= lim ^M= u '(-1), (4.1) 
xj.-i —a; 

and similarly /i(l) = — u'(l) at the right-hand end-point. In short, h G C{D). 



Lemma 4.1. Suppose that f — \A for some open subset A in D. Let u G C(D). Given a G AnD, 
put A n := [a — 77, a + 77] /or 77 > small. Then 
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Proof, (i) Notice that 4n(a-ij, a+??) 7^ for each 77 > 0. Consequently, (fxA v , 1) = m^nA,,) > 
for each 77 > (small) and likewise for (fxA?,, V0- Write 

C/Xy M ) _ "(a) t/Xy 1) + (fXA v : u - u(a)) _ jj(a)(i - u(a)( 2 

(fXA v ,i>) ^{a){fxA v ,l) + {IXA v ^-i>{a)) W V(a)(^(a)+C 2 ) 

where 

^ ._ (fXA„, u-u(a)) 

1 ' (Jxy 1) 

^ . = l/xy ip-i/>(a)) 

2 ' (fXA v , 1) 

Both these last vanish in the limit rj \. and this leads to the identity. 
(ii) From the estimate Q3.4p . for 77 > small, 

and this establishes the limit. □ 
With this preparation in hand we arrive at the crucial (non-)optimality condition. 

Theorem 4.1. Let t G (0, 2/3), m G N and / = xa £ ?7 t . Assume that A satisfies condition 
\3. 61) . Suppose that a, b <E D with a ^ b such that 

(i) h(a) < h(b); 
(ii) a G <9/l; 
(mj 6 G <9A 

TTien i/iere exists f\ G J7 t suc/i i/ia£ </(/) < J{fi). 

Proof. Write / in the form / = 53j=i XA k for some 1 < k < m where Aj — (aj, bj) and 

— 1 < a,\ < b\ < 02 < &2 < • • • < a>k < bk < 1. 
Given 77 > put A (; := [a — 77, a + 77] and B r] := (b — 77, 6 + rf). Set <? := 1 — /. The functions 

V ^ (XA V f, i>) and V ^ (XB C 9, ip) 

are strictly increasing at least for n > small. For e > sufficiently small, there exist unique 
rj > and C > depending upon e such that 

Define 

fe ■= f - fXA„ +9XB ( - 
Then f £ G [/ t (m) for e > small. Now 

J(fe)-J(f) = (fe,Gf e )-(f,Gf) 

= (fe-f,G[f e + f]) 

= (-fXA„ + gXB c , G[2f- fxA n + 9XB ( ] ) 

= 2(gxB ( ,u)-2(fxA ri ,u) + (gxB ( - fXAr,,G[gxB ( - fXA n ]) (4.2) 

where u = G f as usual. Thus, by Lemma 14 . 1 1 (with the help of the Cauchy-Schwarz inequality to 
deal with the cross-terms), 

lime- 1 ! J (fe)-J(f)\ =2( h(b)-h(a)) > 0. 
In particular, there exists e > (small) such that /1 := / e G satisfies J(/i) > ■/(/)• □ 
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5 More on the (non-)optimality condition 

In this section we verify condition (i) in Theorem 14 . 1 1 for some particular configurations / in Ut- 
Lemma 5.1. Let f £ Ut for some t e (0, 2/3) and set u := G f. Then 

(i) u'(-l) + u'(l) = -f D xfdm; 

(ii) u'(-l)-u'(l) = J D fdm. 

Proof, (i) Using tj)' = —x and the integration- by-parts formula, 



x f dm = / ip u dm 
D JD 

= ip'(l)u'(l) - ^'(-l)u'(-l) - J ip"u'dm 
= -u'(l)-u'(-l). 

(ii) This follows from J D u" dm = — j f dm = u'(l) — u'(—l). □ 
Assuming that / ^ 0, define a, b G M. by 

a := inf I x > -1 : (X(-i,x], f) > }, 
1-6 := inf | a: > : (x[i-x,i), f) > 0}; 

so that a 6 [—1, 1) and b £ (—1, 1]. 

Lemma 5.2. For each f S £/ t tyji/i / ^ 0, 

'i(y) = Id i 1 ~ x ) f dm f or each y e I -1 ' a ]>" 

fwj 7i(y) = ^ ( 1 + x ) / dm /or eac/i y £ [6, 1] . 
Proof, (i) Suppose that a = — 1. By (|4.ip and Lemma I57H 
fc(-l) = u'(-l) = (1/2) f (l-x)fdm. 

JD 

Now suppose that a e (— 1, 1). Using integration- by-parts, 

u(y) = / u' dm 
J(-i,v] 

u ip"dm 

(-1,1/1 

= -u'(j/)V'(t/) + u'(-l)V''(-l)+ / u'Vdro 

•/(-i.vl 

= -n'^V'^ + u'C-l) 
as u" = — / = m-a.e. on (—1, a]. For the same reason, 

u'(y) - u'(-l) = I u"dm = {). 
Therefore, 

u(y) = ( i - ) „'(-!) = i±W y ( 1 - a) / dm 
from which the statement is clear. Part (ii) follows in a similar fashion. □ 
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Proposition 5.1. Let f G PUt for some t 6 (0, 2/3). With a, b as in 115.1)) assume that 

(i) -1< a < < -a < b < I; 

(ii) f = 1 m-a.e. on (—a, &). 
17ien h(a) < h(b). 

Proof. From Lemma [ 



h(b) — h(a) — —-. r^- -— f f dm + — - — — / x f dm. 



(l-o)(l + 6) J D J (l-a)(l+b)J D 
Put S := {x e (0, 1) : f(-x) = 1 } as in Lemma El Then 



x f dm = / a; { f(x) — /(— x) } m(dx) + / x f dm 
d Js J(0,1)\S 



xfdm> xfdm=(l/2)(b 2 -a 2 ) 

(0, 1)\S J{-a,b) 



making use of (lij. Thus, 



(2-a + b) / xfdm> (1/2) (2 - a + 6) (6 ^ - or) > 6 -a 2 > (a + b) / /dm 

and hence h(b) — h(a) > 0. □ 
In the next two sections, we show non-optimality of polarised configurations in three broad cases. 

6 Two non-symmetric cases 

Let t G (0, 1/3] and imagine a configuration polarised to the right that charges the left-hand 
interval (— 1, 0) but which is not symmetric under reflection in the origin. We show this is non- 
optimal. 

Lemma 6.1. Let m e N and t G (0, 1/3]. Suppose that f G satisfies the properties 

(*) f - Pf; 

(n) (/, X(-i,o)) > 0. 
(Hi) (f, X(o,i) ) > (/, X(-i,o) )■ 

Then there exists f\ G uj:" 1 " 1 with the property that f\ = P f\ such that J(f) < J(fi). 

Proof. We may assume that / = \A where A satisfies condition (|3.6p . We may then write / in 
the form described at the beginning of the proof of Theorem 14.11 By (ii), a% < 0; and by (i), 
bk > -oi. 

Case (a): —1 < a\ and bk < 1. Then, in fact, — 1 < ai < < — a\ < bk < 1. Put 

k\ := min | j = 1, . . . , k : bj > —a% j. 

Suppose first of all that — aj = bk t ■ Decompose A into its symmetric and non-symmetric parts 
A\ and A 2 as in (|3 . 1 1|) . By (Hi), A 2 ^ 0. Write /i := \A X and f% := xa 2 - By symmetry, 
hf 1 (ai) — hf 1 (—ai). Further, hf 2 (ai) < hf 3 (—a\), this being a consequence of (|2.1[) . Therefore, as 
h = hf t + hf 2 , we obtain h(a\) < h(—ax). The conclusion follows with an application of Theorem 
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If —a\ 7^ fefci then — a\ < bk x and / = 1 on (— 01, b kl ). If k\ = k then h(a{) < h(bk) by Proposition 
15.11 On the other hand, if k\ < k define 

k\ k 

fx ■= ^2 XAj and f 2 := XA 3 ■ 

3 = 1 j=k! + l 

By Proposition 15. 11 hf 1 (ai) < hf 1 (b kl ). It can be seen from the representation in Lemma l5.2l that 
hf 2 is increasing on [—1, a^+i]. In sum, then, h(a\) < h(b kl ). Now apply Theorem 14.11 once more. 

Case (b): —1 < a\ and b k = 1. In this situation, — 1 < a\ < < — a\ < bk = 1. Define k\ as 
before. The case k\ < k may be dealt with in a similar way to case (a) above. So assume that 
k\ = k. As / is polarised to the right, the interval ((— &i)V0, — ai) must sit inside A and so it must 
hold that ak < —a\. In case t G (0, 1/3), it must also hold that < a^. The situation t = 1/3 
and Ofc = forces ( /, X(-i, 0) ) = contradicting (ii). In cither case, therefore, < < — ai and 
fc > 2. 

Consider the function g := \b where B := D \ A. By Lemma T3. 21 g 6 Thus, 

-1 < ai < -an, < < a k < b k = 1, 

and 5 = 1 just to the right of — a k as g is polarised to the left. This situation corresponds to the 
one described at the start of the consideration of this case but for g instead of /. Use the fact 
that h g = 1 — hf . 

Case (c): a\ = —1. Then b k = 1 and a k < — b\ as / = Pf. Apply the arguments in case (a) to 
the function g. □ 

We now take t 6 (0, 1/3) and imagine a configuration that lies entirely in the right-hand interval 
(0, 1) but that has not yet been pushed rightwards to the maximum extent. We again show 
non-optimality. 

Lemma 6.2. Let m £ N and t € (0, 1/3). Suppose that f € U^ m ' satisfies the properties 
(i) (/, X(-i,o})=0; 
(H) if, X(o,co) > 0. 

Then there exists f% e J7 t with the property that fi = Pfi such that J(f) < J(fi)- 

Proof. Again take / = \A where A satisfies condition (|3.6p and suppose / takes the form described 
at the beginning of the proof of Theorem 14.11 By (i), ai > and by (ii), ai < £ t . Therefore 
a\ < b% < 1; for otherwise, if b\ — 1 then 

(/, ip) > (X(oi,i), ^) > (X(? t ,i). V>) = 

Again borrowing the notation of Theorem 14.11 put f\ := XAi and /2 := ^=2^- By Lemma 

ft/ 1 (ai) = -^— / (l-x)dm= b -^-{l-(l/2)(a 1 +b 1 )} 
i_a w(ai,6i) 1-ai 

and 

M & i) = TT7T / (l + x)dm=^f^{l + (l/2)(o 1 + & 1 )}. 

1 + "1 J(ai,bi) l + Ol 

A little algebra yields hf 1 (b\) > hf x (a\). Lemma also indicates that hf 2 is monotone increasing 
on [—1, 02]. Therefore, h(b\) > h{a\). The conclusion now follows with the help of Theorem 

rcn □ 
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7 The symmetric case 

In the last of the three cases, we consider a configuration that is symmetric under reflection in the 
origin. 

Proposition 7.1. Let m G N and t G (0, 2/3). Suppose that f G L^™ 1 ' satisfies the properties 

(i) f = Pf; 

(H) (/, X(-i,o) ) > 0; 
(Hi) if, X(o,i) ) = (/> X(-i,o) )■ 

TTien £/iere exists f\ G L/^*"^ zot£/i the property that fx = P fi such that J(/) < J (fx)- 

Before embarking on the proof of Proposition 17. 11 we require a number of supplementary results. 

Lemma 7.1. Suppose that f = xa £ J7 4 /or some m G N and £ G (0, 2/3). Assume that A 
satisfies condition \3.6\) . Suppose that a G d A H D. Then 

(/xa„, g[/xaJ) . 

')i0 (/XA„, W 

where A v = [a — a + 77] as before. 
Proof. Write 

G[/xaJ(*) = il>(a) (f X A„, 1) + i1>(x) - V(a) } (/XA„ 1) + r(x) 
where r(x) :— ( G(x, ■) — ip(x), fxA n ) for x G -D. Since ^(x) = G(x, x), the estimate f|3 . 2() gives 

|(/XA„r)| <2r 1 (fxA„ I) 2 • 
Forming the inner product we obtain 

(/XA„ G[/ X aJ) = ^(o) (/x^, I) 2 + CfXA,,, i> - V(a)) (/XA„ 1) + (/xa„ r). 
It is clear from this that 



,. (/XA„ , G [fXA v ]) , , 

hm 7T0 = "0(a) 

^0 /XA,, 1 2 



A short step leads to the assertion. □ 

Lemma 7.2. Suppose that f — xa G uf: for some m G N and i G (0, 2/3) and that A satisfies 
condition H3.6]) . Suppose that a, b G D with a 7^ b such that both a G dA and b G dA. Put 
g := 1 — /. Given e > sufficiently small there exist unique r\ > and £ > depending upon e 
such that e = (/ xa,, , VO = (ff Xb c > i>) ■ Then 

lim (/XA,, G[. 9 XB C ]) _ G(a, b) 



eio e 2 ip(a) ip(b) ' 

Proof. Write 

G[gXB ( ](x) = (G(x,-), gxB c ) = G(a, b) (gXB ( , 1) + r(x) 
where r(x) :— (G(x, •) — G(a, b), gxB ( )- For x G A v and y G B^, 
I G(x, y) - G(a, b)\<r, + ( 
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by (|3.2[) . Consequently, 

\(fXA v ,r)\ < (77 + C) (/XA„ 1)(5XB C , 1). 

Now, 

(/XA„, G[gxB ( }) = G(a, b) (fxA n , 1) (ffXs c , 1) + (/XA„, 0. 

from which we derive 

,. (fXA v , G[gxB c }) rv 
urn— —. — r = G(a, 0), 

"4° (/XA„, 1)G?XB C , 1) ^ ^ 

and the conclusion follows straightforwardly. 

□ 

Lemma 7.3. Let b G (0, 1) and a := —6. Lei 77 > small and define £ = £(77) via the relation 

(X[a,a+ V ], ip) = {X{b,b+Q, V0- 
Then £ depends smoothly upon n in a neighbourhood of rj = and 

C = ^+ T ^r ? 2 + 0( ?? 3 ) 

in the limit rj ], . 

Proof. A short computation gives that 

(X[a, a+„] , i>) = (1/2) { (1 - a 2 ) 77 - a ry 2 - (1/3) ry 3 } . 
Define smooth functions /, g : M — > M by 

/fa) := (1 - 6 2 ) 77 + ^ 2 - (1/3) r? 3 and g(Q := (l - 6 2 ) C - K 2 - (1/3) C 3 - 

Now /'(0) = ff'(O) = 1 - 6 2 > 0. In particular, / is strictly increasing in a neighbourhood of 
77 = and g possesses a local smooth inverse h in the neighbourhood of £ = by the inverse 
function theorem. Note that £ is characterised by the relation g(£) = f(rj) for 77 > small. Thus 
£ = {h o /)(??) and depends smoothly upon ?7 for 77 > small. Implicit differentiation yields 
C'(0) = 1 and C"(0) — Trt^- Taylor's theorem with remainder then yields the expansion. □ 

Proof of Proposition ["777} We may suppose that / = \A where A satisfies condition (I3.6[) . Define 
a as in (|5.ip . Assume in the first instance that a e (—1, 0). Put b := —a. Conditions (i)-(iii) 
entail that / is even. In particular, a, b E dA. As in Lemma l?T2l we write 

e = (/XA„, = (sXBc -0) 

for e > small. We aim to show that J(f e ) — J(f) > at least for e > small as in Theorem 14. II 
and we shall borrow notation without comment from its proof. We first claim that (see (|4.2jl ) 

J(f e ) - J(f) 

eiO £ ~ = efo £ ^ 2 ^ 9XBc ' ^ ~ 2 ^ XA7 ' ' ^ + ^ XB< ~ ^ XAv ' G [ 5XBc ~ ^ 64 " ] ^ ^ 

{6/i(6) + tt'(6) + 6(l-6)}. (7.1) 



^(6) 2 

As u e C 1 ' 1 / 2 p) we have that 

(f XAv , u) = «(o) 77 + u'(o) (l/2)ry 2 + O^ 5 / 2 ). 
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Also, with the help of Lemma T7. 31 

[g X B V u) = U (6)C + «'(6)(1/2)C 2 + 0(C 5/2 ) 

= < b ) { v + v 2 } + «'(&) (1/2) v 2 + o( v 5 ' 2 ) 

= < b ) V + { < b ) + \ u '( b ) }*? + 0(v 5/2 )- 

Now u'(a) = —u'(b) because u is even so 

,. (qXb ( , u) - {fxA n , u) ( r],2{9XB ( ,u)-{fxA n ,u) 1 

hm i =hm(-J - = — — - {bh(b) + u (6)} . 

da e 2 eio y e J rj 2 tp(b) 2 

On the other hand, from Lemmas 17.11 and 17.21 

Um (9XB c -fXA„G[gxB ( -fXA,]) 2 G(a, b) = 26(1-6) 

sio e 2 ip(a) tp(a) ip(b) i})(b) 2 

as ip(b) — G(a, b) = b (1 — b). The combination of these identities establishes the claim (|7.1[) . 
We now show that the expression in (|7.ip is positive. From Lemma 15.21 and the even property of 

h ^ = TTbJ D fdm 

and from Lemma 15. 11 

u'(6) = «'(1) = -\ I fdm. 



D 



So we may write 



6 ft(6) +«'(&) + 6 (1-6) = [t^-^} J fdm + b(l-b) 



b 2 (l-b) 
1 + 6 

The conclusion now follows by Theorem 14.11 



> 0. 



The case a = — 1 may be dealt with by applying the above argument (with appropriate modifica- 
tions) to g := 1 — /. □ 

8 The main result 

We are now in a position to prove the main result Theorem 11.31 



Theorem 8.1. Let t e (0, 1/3] and m G N and put A t := (£ t , 1) with £ t as in $1.3)) . Then 
<^ m) = J {f) where f = XA t - 

Proof. By Theorem 13 . 1 1 there exists / <E f7 t such that 0^ — J{f). Now J7 t is closed under 
polarisation. SoP/e t/ t (m) and J(f) < J(Pf) by Theorem ??. We may assume therefore that 
/ = Pf- 

Assume that (/, X(-i,o)) > 0. Since / = Pf it must be the case that (/, X(o, i)) — (/> X(-i,o)) 
(in consequence of Lemma 13. ip . By Lemma 16.11 and Proposition 17.11 there exists /i G with 
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the property that f\ = Pfi and </(/) < J(fi)- This contradicts the optimality of /. We conclude 
that (/, X(_i,o)) = 0. 

lit = 1/3 this compels / = X(o, l) bearing in mind that / £ and / = P f. So let us now take 

t G (0, 1/3). Suppose that (/, X(o,£ t )) > 0- Then the requirements of Lemma l6.2l are satisfied and 

hence there exists fx £ uf: with the property that f\ = P fi such that J{f) < J(fi)- Again this 

contradicts optimality. Hence, (/, X(o,£ t )) = 0- I R iact ' (/> X(-i,ft)) = ®- As / £ L r ( ' m ' we draw 
the conclusion that / = \A t ■ D 

Corollary 8.1. Let t e (0, 1/3]. Then a t = J(f) where f = XA t - 

Proof. Let / £ Ut- By Lindelof's theorem, we may write / in the form / = xa where A = 
Ufcli Ak is a countable union of disjoint open intervals Ak in D. Put /„ :— Ylk=i XA k - By the 
monotone convergence theorem, J(f) = limn^oo J(/„). Note that /„ £ f/ t ■ By Theorem 18.11 
</(/«) < J{XA t ) and J(/) < J(xA t ) on taking limits. □ 

Proof of Theorem ] 1.3[ We only need to deal with the case t £ (1/3, 2/3) in view of Corollary 18. II 
Let / £ Vt for such a i. Then g := 1 — f € Vij-i-t and 

J(/) = 2(i - 1/3) + J(g) < 2(t - 1/3) + J( X A 2/3 _ t ) = J(l - XA 2/3 _ ( ) = J(xaJ. 

This clinches the result in the final case. □ 

9 Application: maximum flux exchange flow 

In this section we prove Theorem 11.41 

Proposition 9.1. It holds that 

(i) 7 A = 2ai_A - (1/3)(1 - A) 2 for A £ (-1, 1); 

(ii) 7 = sup Ae( _ lil) 7A. 

Proof, (i). Fix A G (—1, 1). Let A be an open subset in D. Suppose that u satisfies (|1.4p along 
with the flux-balance condition (u, 1) = 0. Put 

/ = /a, a := 

Then u = Gf. From the flux-balance condition and symmetry of the Green operator, 

= (1, Gf) = ty, /) = -(A + 1)0/;, X a) - (A - 1)W>, X d\a); 
so that A = (tp, iy 1 (ip,XD\A ~ Xa) and (tp, xa) = !)■ Moreover, 

(Xfl\A, u) = (Gxz5\a, /) 

= -( A + i){GXo\A, Xa) - (A - 1)(Gxd\a, Xd\a) 

= 2(^i, 1) _1 | - (V>, Xd\a)(Gxz5\a, Xa) + W>, Xa)(Gxd\ j 4, Xd\/i)} 

= 2(-0, l)" 1 ^, 1)(Gxa, Xa) - Xa) 2 } 

= 2J(xa)-(1/2)(1-A) 2 (V^, 1), 

upon writing Xd\a — Xd — XA on each occurrence in the penultimate line. Now simplify using 
(|3.3[) . These considerations lead to the reformulation (i). The statement in (ii) then follows 
immediately. □ 



-(A + l) on A, 
-(A-l) on D\A, 
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Proof of Theorem\L^ Part (i) follows from Theorem O For t £ (0, 2/3), 

a t = 2(t - 1/3) + a 2 / 3 _t, 
as can be seen from the proof of Theorem 11.31 Therefore, for A £ (—1, 1), 

7_ A = 2ai+A - (1/3)(1 + A) 2 =2ci_ A - (1/3)(1 - A) 2 = 7A . 
We now show that 7a < 70 for each A £ (0, 1). 
Let £ £ (—1, 1) and u :— G f where / := X{£, l)- Then 

u(x)-l 3^-l) 2 (^+ 1 ) if "!<*<£, 

~ I + j(£ + !) 2 1 + 3 - ?(* + X ) 2 if 

A computation leads to 

Also (see (TO)) ). 

^ = (V>, /) = ^) = ^{2-3£ + e 3 }. 
Therefore, 

7A = 2a ¥ -(1/3)(1-A) 2 
= 2J(/)-3^) 2 

12 4^ 4 s 12 s 
=: HO- 

Now ft(Q) = ± an d = and = -(1/2) £(1 - £ 2 ) 2 < for £ e (0, 1). This shows that 

7a < 7o for each A £ (0, 1) as desired. The result follows from this and (i) of the Theorem. □ 
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